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Abstract 

Let Xi denote free identically-distributed random variables. This pa- 
per investigates how the norm of products n„ = XiX2---Xn behaves as n 
approaches infinity. In addition, for positive Xi it studies the asymptotic 
behavior of the norm of Y„ = Xi o X2 o ... o X„, where o denotes the 
symmetric product of two positive operators: Ao B =: A^^^ BA^^^ . 

It is proved that if the expectation of Xi is 1, then the norm of the 
symmetric product Yn is between an^^^ and C2n for certain constant ci 
and C2. That is, the growth in the norm is at most linear. 

For the norm of the usual product Pi„, it is proved that the limit of 
log Norm{Pin) exists and equals log ^JE(XfXi). In other words, the 
growth in the norm of the product is exponential and the rate equals the 
logarithm of the Hilbert-Schmidt norm of operator X. 

Finally, if tt is a cyclic representation of the algebra generated by Xi, 
and if 5 is a cyclic vector, then log Norm{n (II„) ^) — log y/E{XfXi) 
for all n. In other words, the growth in the length of the cyclic vector is 
exponential and the rate coincides with the rate in the growth of the norm 
of the product. 

These results are significantly different from analogous results for com- 
muting random variables and generalize results for random matrices de- 
rived by Kesten and Furstenberg. 



1 Introduction 

Suppose Xi, X2, ... , Xn are identically-distributed free random variables. 
These variables are infinite-dimensional linear operators but the reader may 
find it convenient to think of them as very large random matrices. The first 
question we will address in this paper is how the norm of n„ = XiX2.-.Xn 

behaves. If Xi are all positive, then it is natural to look also at the symmetric 

1/2 1/2 

product operation o defined as follows: Xi o X2 — X^ X2X^ . The benefit is 
that unlike the usual operator product, this operation maps the set of positive 
variables to itself. For this operation we can ask how the norm of symmetric 
products Yn = Xio X20 ...o Xn behaves0 
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Products of random matrices and their asymptotic behavior were originally 
studied by Bellman (1954) One of the decisive steps was made by Furstenberg and Kesten (1960)[ 



who investigated a matrix- valued stationary stochastic process Xi, ... , X„ 
, and proved that the limit of n~^E {\og\\Xi...Xn\\) exists (but might equal 
±00) and that under certain assumptions log |jXi...X„|| converges to this 
limit almost surely. Essentially, the only facts that are used in the proof of 
this result are the ergodic theorem, the norm inequality IIX1X2II < H-'^'ill ||-'^2|| 
and the fact that the unit sphere is compact in finite-dimensional spaces. It is 
the lack of compactness of the unit sphere in the infinite-dimensional space 
that makes generalizations to infinite-dimensional operators non-trivial (see 



Ruelle (1982) for a generahzation in the case of compact operators). More work 



on non-commutative products was done by Furstenberg (1963) Oseledec (1968) 
Kingman (1973)[ and others. The results are often called multiplicative ergodic 



theorems and they find many applications in mathematical physics. For exam- 



ple, see Ruelle (1984) 



In this paper, we study products of free random variables. These variables 
are (non-compact) infinite-dimensional operators which can be thought of as a 
limiting case of large independent random matrices. 

Suppose that Xi are free, identically-distributed, self-adjoint, and positive. 
Suppose also E (Xi) = 1. Then we show that the norm of Yn = Xio X2 o ... o X„ 
grows no faster than a linear function of n. Precisely, we find that 

lim sup n^^ ||y„|| < ci \\Xi\\ . 

n — >oo 

We are also able to show that if Xi is not concentrated at 1, then 
lim inf n"i/2 ||r„|| > C2 > 0. 

n — >QC 

For the usual products n„ — XiX2...Xn we can relax the assumption of 
self-adjointness. So, suppose that Xi are free and identically-distributed but 
not necessarily self-adjoint. Also, we do not require that E {Xi) — 1. Then we 
show that 



n — >oo 



lim log ||n„|| - log Je {X*Xi). (1) 



Another way to describe the behavior of n„ is to look at how the norm of 
a fixed vector ^ changes when we consecutively apply free operators Xi, A"„ 
to it. More precisely, suppose that the action of the algebra of variables Xi on 
a Hilbert space H is described by a cyclic representation tt and that the vector 
^ is cyclic with respect to the expectation E. By definition, this means that 
E (X) — TT (X) ^) for every operator X from a given algebra. Then we show 
that 



ilog||7r(n„)ell =log\/i?(X*A,). (2) 



this convention is not important for the question that we ask. First, it is easy to check that 
Xi o X2 has the same spectral distribution and therefore the same norm as X2 o Xi . Second, 
if Xi, X2, and X3 are free, then the spectral distribution of (Xi 0X2) o X3 is the same as 
the spectral distribution of Xi o [X2 o X3) , and therefore these two products have the same 
norm. In brief, if Xi are free, then the norm of Xi o X2 o ... o Xn does not depend on the 
order in which Xi are multiplied by the operation o. 
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Note that we do not need to take the limit, since the equahty holds for all n. 

The reader may think of cyclic vectors as typical vectors. For example, if 
the representation tt is cyclic and irreducible then cyclic vectors are dense in H. 
In colloquial terms, ^ says that for large n the product n„ cannot increase the 
norm of any given vector ^ by more than [E . And ^ says that for 

every cyclic vector f this growth rate is achieved. 

One more way to capture the intuition of this result is to write 

lim n^"'" log ||n„|| = lim n^"'"log sup ||7r(n„)a;|| 

We have shown that this limit is equal to 

n-ilog||7r(n„)e|| 

where ^ is a cyclic vector. Thus, for large n the product n„ acts uniformly in 
all directions. Its maximal dilation as measured by supn^n^^ ll"" O^n) x\\ has the 
same exponential order of magnitude as the dilation in the direction of a typical 
vector ^. 

It is helpful to compare these results with the case of commutative ran- 
dom variables. Suppose for the moment that Xi are independent commutative 
random variables with positive values. Then, 

lim n-Hog\\Xi...Xj =log||X,||, 

n — »oo 

where the norm of a random variable is the essential supremum norm (i.e., 
||X|| = esssup^^gQ |X (w)!). Indeed, for every e > the measure of the set 
{uj : \Xi (uj) ...Xn > \\Xi\\ ... \\X n\\ ~ e\ is positive. Therefore ||Xi...X„|| = 
IjXill". Note that log (X*X^ < log||Xi|] and therefore the norm of free 
products grows more slowly than we would expect from the classical case. 

Another interesting comparison is that with results about products of ran- 
dom matrices. Let Xi be i.i.d. random k x k matrices. Then under suit- 
able conditions, hm„^oo log || exists almost surely. Let us de- 



note this hmit as A. Furstenberg (1963) developed a general formula for A, and 
Cohen and Newman (1984) derived explicit results in the case when entries of 
Xi have a joint Gaussian distribution. In particular, if all entries of Xi are inde- 
pendent and have the distribution TV (O, si) then A = (1/2) {log (s^) -I- log 2 -I- {k/2)} 
where ip is the digamma function (-0 (a;) = d log T (x) / dx) . If the size of the ma- 
trices grows (fc oo) then A ^ (1/2) log (fcs^) . To compare this with our 
results, note that if ks^ —>■ s^, then the sequence of random matrices approxi- 
mates a free random variable Xi with the spectral distribution that is uniform 



inside the circle of radius s. For this free variable, E = and 



our 



theorem shows that lim^ 



Xi...Xr, 



— logs. This limit agrees with the 



result for random matrices. Thus, our result can be seen as a limiting form of 
results for random matrices. 
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The results regarding are also interesting. We can associate with Xi 

and Yn probability measures /i^ and /iy^, which are called the spectral proba- 
bility measures of Xi and Y„, respectively. Then the measure fiy^ determined 
only by n and the measure fj^x and is called the n-time free multiplicative con- 
volution of fix with itself: 

/iy^ = II X ... H fix- 



The norm \\Yn\\ is easy to interpret in terms of the distribution fiy^- Indeed, 
it is the smallest number t such that the support of /iy is inside the interval 
[0, t] . Therefore, the growth in ||K„|| measures the growth in the support of the 
spectral probability measure if the measure is convolved with itself using the 
operation of the free multiplicative convolution. 

In the case of classical multiplicative convolutions of probability measures, 
the support grows exponentially, so that if fix is supported on [0, Lx] , then the 
measure iix-^ x„ supported on [0, (ix)"] • What we have found in the case 
of free multiplicative convolutions is that if we fix EX^ = 1, then the support 
of the fiY„ grows no faster than a linear function of n, i.e., the support of fiy^ 
is inside the interval [0, cnLx] with an absolute constant c. 

As was pointed out in the literature, a similar phenomenon occurs for sums 
of free random variables. The support of measures obtained by free additive con- 
volutions grows much more slowly than in the case of classical additive convolu- 



tions. This effect was called superconvergence by Bercovici and Voiculescu (1995) 
Our finding about \\Yn\\ can be considered as a superconvergence for free mul- 
tiplicative convolutions. 

The rest of the paper is organized as follows. Section 2 formulates the results. 
Section 3 contains the necessary technical background from free probability 
theory. Sections 4, 5, and 6 prove the results. And Section 7 concludes. 



2 Results 

A non- commutative probability space (.4, E) is a unital C*-algebra A and a 
positive linear functional such that E {I) = 1 . We will assume that the 
functional is tracial, i.e., E {AB) = E {BA) for any two operators A and B 
from algebra A. The elements of algebra A are called random variables and the 
functional E is called the expectation. The numbers E {X'^^ are called moments 
of the random variable X. 

A prototypical example of a non-commutative probability space is a group al- 
gebra. That is, for a countable group G we consider the Hilbert space (G, v) , 
where v is a counting measure, and consider the left action of G on (G, v): 
if / € (G, and a,b G G, then [af] (b) = f{ab). The elements of the group 
algebra A are finite sums X^aeG ^<^^ '^^^ extend by linearity the action 

of the group G on (G, v) to the action of the algebra A on (G, v). We can 
additionally complete the resulting operator algebra in an appropriate topology. 
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The expectation of an element '^XaO, is defined as Xe, where e is the identity 
of the group. 

Another important example is the algebra of random matrices. The expec- 
tation of an element X in this algebra is defined as E {X) — £ (^N^^ tr (X)) , 
where £ is the expectation with respect to underlying randomness and iV is the 
dimension of the random matrix. For more details about these examples the 



reader may consult Hiai and Petz (2000) 



The concept of freeness substitutes for the concept of independence. Con- 
sider sub-algebras. Ai^..., An be given. Let at are elements of these sub-algebras 
such that tti £ Ak(i)- 

Definition 1 The algebras Ai^..., An (and their elements) are free, ij E {ai...am) 
0, provided that E (oj) = 0, k(i) =/= k {i + I) for every i < m, and k (to) ^ k (l) . 

Consider the group algebra for a free group with at least two generators. 
Then the operators corresponding to generators are free in the sense of the 
previous definition. For the algebra of large random matrices, Voiculescu proved 
the asymptotic freeness of two classically independent Gaussian matrices, where 
asymptotic means that the property in the previous definition is approached as 



the dimension of matrices N ^ oo (see Voiculescu (1991) I 



It turns out that many concepts of classical probability theory can be trans- 
ferred to the case of free random variables. For example, for a self-adjoint 
variable we can define its distribution function. Indeed, if A is a self-adjoint 
operator then by the spectral decomposition theorem it can be written as 

/oo 
XP (dX) , 
-oo 

where P is a positive, projector- valued measure, i.e., a mapping that sends sets 
of the real axis to orthogonal projectors. This allows definition of the spectral 
measure of A, /x^, which is a measure with the following distribution fimction: 



TA{t)^Ef f P (dX) 



We can calculate the expectation of any summable function of a self-adjoint 
variable A by using its spectral measure: 

/oo 
fWdfiA (A), 
-oo 

Let Xi, X2, ... , Xn be free identically-distributed positive random variables. 
Consider n„ = XiX2.-.Xn and y„ = XioX20...oXn (by convention we multiply 
on the left, so that, for example, X10X20X30X4 = Xio(X2 o (X3 o X4)) ). We 
will see later that these variables have the same moments: E (Iln)*^ = E (1^)*^ . 
As a first step let us record some simple results about the expectation and 
variance of Yn and n„ . We define variance of a random variable A as 

(A) ^■.E{A*A) - \E{A)\\ 
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Proposition 1 Suppose that Xi are self-adjoint and E (Xi) = 1. Then E {Jin) = 

E (r„) = 1 and ^2 (n„) = (y^J ^ _ 

Note that the Unear growth in the variance of n„ — Xi...Xn is in contrast 
with the classical case, where only the variance of log {Xi...Xn) grows linearly. 
We will prove this Proposition later when we have more technical tools available. 
Before that we are going to formulate the main results. 

Let 11 denote the usual operator norm of operator A. 

Theorem 1 Suppose that Xi, X„ are identically- distributed positive self- 
adjoint free variables. Suppose also that E (Xj) = 1. Then 

(1) there exists such a constant, c, that ||1^|| < c||Xi|| n; 
and 

(2) ||r„|l >a(XOv^. 

For the next theorem define 

Theorem 2 Suppose that Xi, X„ are free identically- distributed variables 
(not necessarily self-adjoint). Then 

(1) there exists such a constant, c, that ||n„|| < c ||Xj|| [£■ (X*Xi)]^"~^^^^ ; 
and 

(2) ||n„|| >7VV/4 [E{X*Xi)f\ 

Corollary 1 Suppose that Xi, X„ are free identically- distributed variables 
(not necessarily self-adjoint). Then 

lim n-i log ||n„|| = log jE{X*Xi) 

n— >oo V 

Next, suppose; that the algebra A acts on an (infinitely-dimensional) Hilbert 
space H. In other words, let tt be a representation of A. We call representation 
TT cyclic if there exists such a vector ^ G H that E (X) = (^, n (X) ^) for all 
operators X G A. The vectors with this property are also called cyclic. 

Theorem 3 Suppose n is a cyclic representation of A, ^ is its cyclic vector, 
and Xi, X„ are free identically- distributed variables from A. Then 

log IItt (n„) Cll = log ^JE{X*Xi) 
Corollary 2 If tt and ^ are cyclic then 

log ||n„|| ~ log IItt (n„) ^11 ~ nlog ||7r (Xi) ^|| 

as n ^ oo. 
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3 Preliminaries 

The Cauchy transform of a bounded random variable A is defined as follows: 

, 1 \ 1 M'=^ 



A J z ^ z'^+i 

This power series is convergent for \z\ > . Let us also define the ip-function 
of A: 

^ ^ fe=i 

The i/^-function is convergent for \z\ < \\A\\^^ and it is related to the Cauchy 
transform by the following equality: 

G„ (z) = (z-i) + 1] . 

If A is bounded and E (A) ^ 0, then for z in a sufficiently small neighborhood 
of 0, the inverse of V'a (z) is defined, which we denote as ip^^ (z) . Then the 5- 
transform is defined as 

Sa (z) = (l + ^) V-a' (^) (3) 

Let us write out several first terms in the power expansions for ijj (z) , 
tlj~^ (z) , and S (z) . Suppose for simphcity that E {A) = 1 and let E (A*^) — nik- 
Then, 

i/j (z) — z + m2Z^ + msz^ + 

(z) = z — m2Z^ — (m3 — 2TO2) z^ + 
S'(z) = 1 + (1 — ?7i2)z + (2m2 — m2 — ms) z^ + ... 

The main theorem regarding the multiplication of free random variables was 
proved by [Voiculescu (1987) Later the proof was significantly simplified by 



Haagerup (1997) 



Theorem 4 (Voiculescu) Suppose X and Y are bounded free random vari- 
ables. Suppose also that E (X) 7^ and E {Y) ^ 0. Then 

SxY {z) = Sx {z) Sy (z) . 

In particular, this theorem implies that Su^ = Sy^^ = {SxY^ , where Sx 
denotes the S'-transform of any of Xi . Now it is easy to prove Proposition [TJ 
Indeed, let us denote iS'n„ as Sn- Then, using the power expansions we can write: 

Sn{z) = 1+ fl-m^"^)z + ... 



'2 

= (5'x)" = 1 + ?^ (1 - "12) 2 + 

where m.j"'' —'■ £- (II„)^ and TO2 —: E [Xi)^ . Then, using power expansion in 
([3]), we conclude that i?(n„) = 1. Next, by definition, {Xi) = m2 — 1 and 
(7^ (n„) = 7712"^ — 1. Therefore, we can conclude that (n„) = na'^ {X) . QED. 
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4 Proof of Theorem [T] 



Throughout this section we assume that Xi are self-adjoint, E (Xi) — 1, and 
the support of the spectral distribution of Xi belongs to [0, L] . 

Let us first go in a simpler direction and derive a lower bound on ||i^Ti|| • 
That is, we are going to prove claim (2) of the theorem. From Proposition 
m we know that E (Yn) = 1 and (Yn) = na^ (Xi) . It is clear that for ev- 
ery positive random variable A, it is true that E {A^) < \\A\\ and therefore 



ll^ll > yo-M^) + [E{A)f. Applying this to r„, we get ||r„|| > Vncr^ + l. In 
particular, ||F„|| > a^/n, so (2) is proved. 

Now let us prove claim (1). By Theorem HI Sn{z) = {Sx {z))". The idea 
of the proof is to investigate how l^x {z)\" behaves for small z. It turns out 
that if z is of the order of then \Sx (z)\" > c where c is a constant that 
does not depend on n. We will show that this fact implies that (z) (i.e., the 
V'-function for Yn) has the convergent power series in the area \z\ < (cn)~^ and 
that therefore the Cauchy transform of Yn has the convergent power series in 
\z\ > cn. This fact and the Perron-Stieltjes inversion formula imply that the 
support of the distribution of F„ is inside [—cn, cn] . 

In the proof we need the result about functional inversions formulated below. 
By a function holomorphic in a domain, D, we mean a function which is bounded 
and differentiable in D. 

Lemma 1 (Lagrange's inversion formula) Suppose f is a function of a com- 
plex variable, which is holomorphic in a neighborhood of zq = and has the 
Taylor expansion 

oo 

f{z) ^ z + J2'^kz'' + 

k=2 

converging for all sufficiently small z. Then the functional inverse of f (z) is 
well defined in a neighborhood of and the Taylor series of the inverse is given 
by the following formula: 



k=2 



where j is a circle around 0, in which f has only one zero.. 



For the proof see Theorems II. 3. 2 and II. 3. 3 in Markushevich (1977) 



Section 7.32 in |Whittaker and Watson (1927) 
Lemma 2 E {X'') < L^-^ . 
Proof: 

E{X'')= A'^d/i^ (A) < L'^-M Ad/^^ (A) = 
where dfix denotes the spectral distribution of the variable X. QED. 



Lemma 3 The function ipx (z) is has only one zero in \z I < (4L) , and 
\z\ = (4L)"^ , then IV-^ (z)| > (6i)"^ . 

Proof: If |2| < (4L)"^ then 

OO 

\i>x{z)-z\ < \z\J2e{x'')\z\''~' 

1 \z\ 



k^2 

OO 



< NE 



k=l 



4fc 3 



Therefore, by Rouche's theorem, i/'x {z) has only one zero in this area. 
If \z\ = (4L)"^ , then 

OC 

IV-xWI > \z\-Y.E{x'')\zf 

k=2 



> Nil 

4L 



1 - 



4fe 

/c=l / 

1\ 1 



3/ 6L 



QED. 



By Lemma [U we can expand the functional inverse of ■(/'x i^) ^ follows: 

OO 

where 



k=2 



Ck 



2'Kik 



dz 



1 [^x (2)] 



Lemma 4 // |u| < (72Ln) ^ , t/ien 



< 



7n 



Proof: Using the previous lemma we can estimate Ck'- 



Then 



< 



k=2 



fc=l 



12n 



3 1 



2 12?^ - 1 



3 12?i 1 1 
2 12n - 1 12^ - 7^' 
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provided that |u| < (72Ln) ^ . QED. 
Lemma 5 If\u\ < {72Ln)~^ , then 

Proof: Recall that Sx {u) = (1 + u) V'x^ («) /u- Then we can write: 



\^-Sx{u)\ 



< \u\ + \l + u 



(I+.)ffi^-l 



Then the previous lemma implies that for |m| < (72Ln) ^ and n > 2, we have 
the estimate: 

|l-^-(«)l^72k + 
Note that L > 1 because EX — 1. Therefore, 



1 



72Ln 



1 

7n' 



QED. 



Lemma 6 For all positive integer n if \u\ < (72Ln) ^ , then 

Proof: Let us first prove the upper bound on \Sx . The previous lemma 
implies that 

Now let us prove the lower bound. The previous lemma implies that 



In an equivalent form, 



n\og\Sx {u)\ > nlog 1 - 



6n 

Recall the following elementary inequality: If a; € [0, 1 — e~^] , then 

log {l-x)> -2x. 

Let a; = 1/ (6n). Then 

6nJ 3n 



(4) 
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Substituting this in ([3]), we get 



nlog\Sxiu)\ > 



3' 



or 



\Sx{ur>e-'/'. 



QED. 

By Theorem |4l S'„(u) —: [Sx (u)]" is the S'-transform of the variable Yn. 
The corresponding inverse -0-function is (u) = uSn (u) / (1 + u) . 
First, we estimate Sn (u) — 1. 

Lemma 7 // |u| < (72Ln)^^ , then 



\Sn{u)-l\<-. 



Proof: Write 

\sxiur-i\ < \sx{u)-i\(\sx{u)r' + \sx{u)r' + ... + i) 



6n 5 



QED. 



Lemma 8 The function tp^^ (u) has only one zero in \u\ — (72Ln) ^ and if 
\u\ — {72Ln) , then 

U-i (u)\ > 

1^" ^ - 102Ln 



Proof: Recall that by definition in ([3]), V'n^ (") — uSn (u) / {1 + u) . There- 



fore, 



Sn{u)-il+u) 



1 + u 



and by Lemma [7] we have the following estimate: 



Sn (u) - (1 + U) 



1 + U 



< -^\Sniu)-l\+ 



1-H 

72 1 1 1 

< \ < -. 

- 71 5 71-4 



Therefore, by Rouche's theorem, ip„^ (u) has only one zero in |u| < (72Ln) ^ . 
Next, note that ip~^ (u) = uSn (u) / {1 + u) and if \u\ = {72Ln)^^ , then 



1 + u 

Using Lemma m we get: 



> 



1 



72Ln' 



/ 1 + 



1 



> 



1 



72Ln ~ 73Ln 
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QED. 

Now we again apply Lemma [T] and obtain the following formula: 



k=2 L 



If du 



1 [i^n^u)] _ 



Z 



(5) 



where we can take the circle \u\ — (72Ln) ^ as 7. 

Lemma 9 The radius of convergence of series ^ is at least (102in) ^ . 

Proof: By the previous lemma, the coefScient before z'^ can be estimated 
as follows: 

This implies that series ([S]) converges at least for \z\ < (102Ln)^^ . QED. 

Lemma 10 The support of the spectral distribution of Yn = Xi o X2 o ... o X„ 
belongs to the interval [— 102Ln, \Q2Ln] . 

Proof: The variable K„ is self-adjoint and has a well-defined spectral mea- 
sure, /i„ {dx) , supported on the real axis. We can infer the Cauchy transform 
of this measure from (z): 

Gn {Z) = + 1] . 

Using Lemma ini we can conclude that the power series for Gn {z) around z = 
00 converges in the area \z\ > 102Ln. The coefficients of this series are real. 
Therefore, using the Perron-Stieltjes formula we conclude that /i„ (dx) is zero 
outside of the interval [—102Ln, 102Ln] . QED. 
Lemma [10] implies the statement of Theorem [TJ 



5 Proof of Theorem [2] 

The norm of the operator n„ coincides with the square root of the norm of 
the operator n,*n„. Therefore, all we need to do is to estimate the norm of the 
self-adjoint operator n*n„. 

Lemma 11 For every bounded operator X (z A, products X*X and XX* have 
the same spectral distribution. 

Proof: Since E is tracial, E{X*X)'' = E{XX*)'". Therefore, X*X and 
XX* have the same sequence of moments and, therefore, the same distribution. 
QED. 

If two variables A and B have the same sequence of moments, we say that 
they are equivalent and write A ~ _B. In particular, two self-adjoint bounded 
variables have the same spectral distribution if and only if they are equivalent. 
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Lemma 12 Let A, B, and C he three hounded operators from a non- commutative 
prohahility space A. If A B, A is free from C , and B is free from C, then 
A + C B + C, AC BC, and CA - CB. 

Proof: Since A and C are free, the moments of A + C can be computed 
from the moments of A and C. The computation is exactly the same for B + C, 
since B and C are also free. In addition, we know that A and B have the 
same moments. Consequently, A + C has the same moments as B + C, i.e., 
A + C ^ B + C. The other equivalences are obtained similarly. QED. 

Lemma 13 If A ^ B, then Sa {z) — Sb {z). In words, if two variahles are 
equivalent, then they have the same S -transform. 

Proof: From the definition of the function, it is clear that ii A ^ B, then 
i^) ~ i-^)- This implies that tp^^ (z) = ■0^^ (z) and therefore Sa {z) = 
SBiz). QED. 

For example, Sx'Xi {z) does not depend on i and wc will denote this func- 
tion as Sx'X [z) . 

Lemma 14 If Xi, Xn are free, then 

n^n„ ^ x*Xn---XiXi 

and if Xi, Xn are in addition identically distrihuted, then 

<S'n'n„ = 5'n„n* = {Sx'x)"' 

Proof: We will use induction. For n = 1, we have n^IIi = X^Xi. Therefore 
'S'n*ni = Sx-x- Suppose that the statement is proved for n—1. Then 

n*n„ = x*-^...xiXi...Xn 

~ XnX*-^X*-^_-^...XlXi...Xn-l, 

where the equivalence holds because E is tracial and it is easy to check that the 
products have the same moments. Therefore, 

n:n„ ^ (x„x:)n:_in„_i 

by Lemmas [11] and 1121 Then the inductive hypothesis implies that 

n*n„ ~ XnXn...xiXi. 

Using Lemma fT3l and Theorem |4j we write: 

'5'n*n„ = {Sx-'xY ■ 

Since n*n„ - n„n;, therefore, S^n-n,. S'n„n- = [Sx-xT ■ QED. 

We have managed to represent Syi*^yi„ as [Sx'xY^ and therefore all the 
arguments of the previous section are applicable, except that we are interested in 
[Sx'xT" rather than in (S'x)" • In particular, we can conclude that the following 
lemma holds: 
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Lemma 15 Define 



7 = cr 



E {X*X, 



Then 

(1) ||n;n„|| < 102\\X£ nE(X*X^f-\ and 

(2) \\wM\>i^E{x:x,T . 

Proof: Let us introduce variables Ri = s^^Xi, where — E {X*X) . Then 
\\R*Ri\[ = (\\X,\\/sf and E{RlRi) = 1. Let n„ = Ri...R„. Then n;n„ = 
s2"n;n„ and the^S'-transform of n;n„ is {Sr,*bT ■ 

Note that n*n„ has the same S'-transform and therefore the same spec- 
tral distribution as (R^Ri) o ... o (R^R^) . Using Thcorcm[Tl we conclude that 
n^fin < 102 (||X,|| /s)^n. It follows that ||n;n„|| < 102||X,f s2"-2„_ In 
addition, Theorem [1] implies that 

n;n„ > V^2<j{r*r,) 



Consequently, 



2n 



QED. 

From Lemma [15] we conclude that 



|n„|| < u\\x,\\V^[E{x*x,)] 



(n-l)/2 



and 



l|n„||>7'/V/4 [E{x:x,)f' 

This completes the proof of Theorem [2l 



6 Proof of Theorem [3] 

By definition of the cyclic vector, we have: 

||7r(njci|2 = (^(n„)^,7r(n„)0 
= (e,7r(n:n„)0 
= i?(n:n„). 

Using Lemma [131 we continue this as follows: 

E{niii„) - E{x:Xr,...xiXi) 

= [E{X*X)Y\ 

Consequently, 



n-Hog\\naC\\^^logE{X*X) 



QED. 
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7 Concluding Remarks 



We have investigated how the norms of n„ = Xi...Xn and Yn = Xi o ... o Xn 
grow as n ^ oo. For ||n,i|| , wc have shown that hm„^oo log l|n„|| exists 
and equals log E {X*Xi). For ||F„|| , we have proved that the growth rate of 
is somewhere between -^/n and n. There remains the question of whether 
lim„_^oc fi^"' ll^n II exists for some s. 

Another interesting question, wliich is not resolved in this paper, is how 
the spectral radius of n„ grows. Indeed, for y„, the norm coincides with the 
spectral radius. But for n„, the norm and the spectral radius are different 
because n„ is not self-adjoint. 
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